Crystal structure and electronic states of tripotassium picene 
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The crystal structure of potassium doped picene with an exact stoichiometry (K3C22H14, K^picene from here 
onwards) has been theoretically determined within Density Functional Theory allowing complete variational 
freedom of the crystal structure parameters and the molecular atomic positions. A modified herringbone lattice 
is obtained in which potassium atoms are intercalated between two paired picene molecules displaying the two 
possible orientations in the crystal. Along the c-axis, organic molecules alternate with chains formed by three 
potassium atoms. The electronic structure of the doped material resembles pristine picene, except that now the 
bottom of the conduction band is occupied by six electrons coming from the ionized K atoms (six per unit cell). 
Wavefunctions remain based mainly on picene molecular orbitals getting their dispersion from intralayer edge 
to face CH/7T bonding, while eigenenergies have been modified by the change in the electrostatic potential. The 
small dispersion along the c axis is assigned to small H-H overlap. From the calculated electronic density of 
states we expect metallic behavior for potassium doped picene. 
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I. INTRODUCTION 



Recently, Mitsuhashi et al. (Ref. 1) have shown the exis- 
tence of superconductivity at relatively high temperatures for 
alkali-metal-doped picene. In particular, K3 ^picene shows a 
critical temperature of T c = 18 K. This work opens the door 
to a new class of materials that promise improved T c values. 
While the structure of the organic picene crystal is well known 
since 1985 2 , the best possible structural determination using 
synchroton radiation on K doped samples has been indecisive 
about the actual positions of interstitial atoms^. No further 
advance in the understanding of the superconductivity is pos- 
sible until a successful detailed structural determination ex- 
ists. Once the atomic structure is known, both the electronic 
states and the phonons can be quantitatively studied and the 
usual ideas to understand superconductivity followed. A pre- 
liminary theoretical determination of the crystal structure of 
pristine and potassium doped picene (K^picene) has been re- 
ported by Kosugi et al. 4 by a first-principles calculation which 
adheres strictly to experimentally determined crystal parame- 
ters for the unit cell vectors length (see also comment in Ref. 
5 ). Kosugi et al. find that potassium doped picene molecules 
show an edge-to face intralayer ordering but stay almost per- 
pendicular to the ab plane, losing the conventional herring- 
bone interlayer stacking (see Fig. 4a of the referred paper). 
In our work, to circumvent the difficulties found in the experi- 
mental determination of the structure of K^.gpicene (small and 
disordered samples, possible existence of more than one struc- 
tural phase, etcpJ) we have adopted a more general strategy 



to the structural search of the crystalline structure of highly 
doped picene; we rely on the best available ab initio meth- 
ods based in Density Functional Theory (DFT) to get a stable 
structure that provides at least a local minimum of the free 
energy of a precise stoichiometry. We study the case of three 
potassium atoms per picene molecule because it is the one 
showing the higher critical temperature. Therefore, in our cal- 
culation atomic positions and parameters defining the unit cell 
(including the symmetry group) have been fully optimized to 
reach a total energy minimum configuration subject to mini- 
mum forces and stresses. The resulting crystalline structure 
shows enough internal consistency and symmetry to be the- 
oretically appealing. We remark that the calculated struc- 
ture differs significatively from two different structures pub- 
lished for a smaller K doping of a similar material (pentacene 
crystaljMl 

The rest of the paper is organized as follows: Section II 
gives details for the computational procedure followed in this 
work. Section III presents our optimized crystalline structure 
and the corresponding electronic band structure and density of 
states. Finally, Section IV summarizes our conclusions. 

II. NUMERICAL PROCEDURE 

A full experimental determination of the doped picene crys- 
tal structure based in X-rays diffraction data has not been pos- 
sible so far because of the large background and the small 
number of useful measured peaks. However, a basic Le Bail 
determination of the structure has been performed to obtain 
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FIG. 1: (Color online) Upper panel shows the optimized atomic 
structure of K3picene as obtained with CASTEP for a completely 
free lattice structure. The P2i spatial symmetry of pristine picene is 
conserved. Notice that although only one atomic layer is displayed, 
the positions of the next-layer picene molecules are on top of potas- 
sium atoms. Lower panel shows an alternative structure obtained 
starting the geometrical optimization from a laminar arrangement of 
picene molecules intercalated by K atoms in a global configuration 
based on intercalated graphite. 



approximate values for the unit cell parameters 3 , yielding val- 
ues quite similar to the ones derived for pristine picene. There- 
fore, we propose to complement the presently limited experi- 
mental information with a pure numerical procedure based in 
Density Functional Theory to define realistic candidates for 
the actual crystalline structure. One important advantage is 
that Poly cyclic Aromatic Hydrocarbons (PAH) show a robust 
geometry even in the presence of alkali metals and translate 
and rotate almost rigidly during the total energy minimiza- 
tion. However, the position of the dopant atoms is a complex 



problem not fully understood, as shown by previous theoret- 
ical studies concluding different alternative models^! In or- 
der to get a plausible crystalline structure for K 3 picene, we 
have minimized the free energy per unit cell using a DFT 
scheme at the Local Density Approximation (LDA) level. 
Both atomic positions and cell parameters (a, 6, c, a, (3, 7) are 
free to change to minimize forces on atoms and stresses on 
the unit cell get. Some more technical parameters of the com- 
putation follow. We use norm-conserving pseudopotentials 9 , 
an energy cutoff of 660 eV, and Monkhrost-PacW^ meshes of 
4x6x3 or 5x5x3 depending on the unit cell being under 
consideration. Actual calculations have been performed with 
the CASTEP program allowing for spin polarization of differ- 
ent bands of the doped crystal 11 12 . To describe the exchange 
and correlation potential we use Kohn and Sham's local den- 
sity approach as computed by Ceperley and Alde J 13 ! 14 ' . Maxi- 
mum forces and stresses on the system are converged to target 
thresholds of 0. 1 eV/A and 0. 1 GPa. 

The whole procedure has been applied to pure picene as 
an overall check. The optimization was started using the 
complete structural determination provided by the Cambridge 
Crystallographic Data Centre (CSD_CIF_ZZZOYC01). Cell 
parameters (a, 6, c, f3) change from (8.480 A , 6.154 A , 13.515 
A , 90.46°) to (8.373 A , 6.043 A , 13.390 A , 90.18°) improv- 
ing the cell total energy in merely 0.12 eV. The unit cell vol- 
ume shrinks about 4%, but the change in the molecule struc- 
ture is imperceptible. Therefore, the test is quite satisfactory 
and we assume that the differences between experimental and 
computational values provide an estimate for the errors in- 
volved in our structural search. 



III. RESULTS AND DISCUSSION 

A. Crystal Structure 

Previous work on K-doped pentacene^ and on Kapicene^ 
has proven that when large enough amounts of alkali atoms 
are introduced in a herringbone organic structure, the pre- 
ferred positions for the dopants is intralayer, i.e., in between 
planes defined by PAH molecules. The other possibility, K 
atoms in the planes between PAH layers is energetically less 
favorable but it has also been considered for Kipentacene 7 . 
Since we are mainly interested in the case of very high doping 
(K 3 picene) because it shows promising results as a new su- 
perconducting material, we have assumed from the beginning 
that K occupies sites between picene molecules. As a start- 
ing crystal structure, we use pristine crystalline picene with 
K atoms in the centroids of next-nearest-neighbors parallel 
picene molecules 15 . Subsequently, the system is allowed to 
freely evolve under atomic forces obtained from a series of 
fully converged DFT calculations to drive the molecules to 
equilibrium positions. We have performed three sets of dif- 
ferent calculations: (i) we froze the unit cell parameters to 
the values obtained from a Le Bail X-ray analysis, (ii) values 
similar to the ones deduced from a Le Bail X-ray analysis, but 
relaxed by less than 10% in order to diminish the structural 
stress, and (iii) fully optimized values for the unit cell lengths 
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FIG. 2: (Color online) Calculated powder X-ray diffraction pattern of 
the optimised theoretical structure of pristine picene (upper panel), K 
doped picene in the optimized herringbone structure (medium panel), 
and K doped picene in the alternative laminar structure (lower panel). 
Intensity is given in arbitrary units while diffraction angle is mea- 
sured in degrees. 



and angles. The optimum configuration for the first case al- 
lows us to minimize the maximum component of forces on 
any atom below F max = 0.03 eV/A, but the maximum com- 
ponent of the tensor stress is S max = 5 GPa, and the total 
energy is above the global minimum found later by 1.42 eV. 
The second case allows us to decrease S ma x = 2 GPa, but the 




FIG. 3: (Color online) Brillouin Zone obtained for K doped picene 
showing the crystalline unit cell and molecules therein. The TFQZr 
path followed to calculate the band structure is shown in red. 



total energy is still above the best configuration by 0.97 eV. Fi- 
nally, the fully optimized structure (a = 7.359A b = 7.36lA 
c = 14.018A a = 90° /3 = 105.71° and, 7 = 90° ) allows us 
to decrease the maximum force and stress below 0. 1 eV/A and 
0.1 GPa respectively. The LDA cell volume increases from 
669 A 3 for clean picene, to 730 A 3 after doping with potas- 
sium. Even after a full optimization the system belongs to the 
space group P2i and remains monoclinic within a ±0.01 A 
error in the atomic positions. Also a = b at the precision of 
our computation. A more detailed analysis shows that picene 
molecules acquire a small regular curvature that contrasts with 
the structure of the free molecule which is completely pla- 
nar. Fig. [ij shows how K atoms are intercalated between the 
picene molecules that form a layer. The high-symmetry of the 
arrangement is clear. Precise structural information is given 
in a supplementary «/-like 

To understand the effect of different starting atomic config- 
urations on the final equilibrium structure, we have looked at 
a laminar structure resembling intercalated graphite. The unit 
cell is made with two picenes and six potassium atoms too, 
but now the picene molecules are parallel. Potassium goes 
to centroids between corresponding hexagonal rings. Start- 
ing from this configuration, DFT forces lead to a locally sta- 
ble structure in which picene molecules form a global plane 
whereas intercalated potassium completes a structure resem- 
bling the well known KCs- The total energy of this structure 
is more than two eV per cell higher than the previous herring- 
bone one. Consequently, from a strict energetic point of view, 
it is clearly inferior. Fig. [TJ) shematically shows this laminar 
structure. Again, precise structural information is given in a 
supplementary ci/-like 

The ultimate justification of a crystal structure should come 
from experiment, typically from X-ray diffraction. Mitsuhashi 
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FIG. 4: (Color online) Density of states and bands around the Fermi level for pristine picene as described by the LDA DFT scheme use in this 
work (upper panels) and the corresponding one after doping with potassium and relaxing the herringbone original structure (lower panels). 
The top of the valence band of both compounds has been aligned to facilitate the comparison. Solid (red) lines indicate the Fermi level position 
for the metallic system. 



et alP provide the X-ray diffraction pattern of a K2.9picene 
sample obtained with synchroton radiation of A = 0.99937A 
in the supplementary information of their paper. We have used 
the mercury program 17 to calculate a powder pattern for any 
well defined crystal structure. Our results are shown in Fig. 
[2] First, we have checked that both the original picene crystal 
and the one refined by DFT at the LDA level provide spectra 
similar to the one given in Ref. 1 Second, we have used the 
same program to get the pattern for our predicted K^picene 
structure (see Fig. [2]). It can be observed that the small angle 
part of the spectrum of doped picene resembles the one ob- 
tained for pristine picene in the herringbone structure. This is 
easily explained since the arrangement of picene molecules is 
still close to the original one. Until here there is some agree- 
ment with experiment since Mitsuhashi et al. comment that 
X-ray diffraction pattern of K 2 .gpicene shows a good corre- 
spondence with that of pristine picene (Figs. S7a and S7b 



of Ref. If- Nevertheless our well ordered array of K atoms 
produces important peaks at larger angles that have not been 
experimentally observed. We believe that K disorder, phase 
mixing, etc. prevent a better experimental determination of 
dopant positions and demand more research in the future. 

On the other hand, calculated powder pattern spectra allow 
us to discard alternative laminar stuctures because they sys- 
tematically show the extintion of the (0,0,1) low angle peak. 
It is easy to check that the spectrum of potassium intercalated 
graphite (KCg) shows a similar property. 



B. Electronic Structure 

The motivation behind the structural work is to better un- 
derstand the electronic properties of the material. We com- 
pute the electronic density of states and the band structure 
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FIG. 5: (Color online) Detail of the conduction bands occupied by 
K doping. Upper panel shows the bands of pristine picene and lower 
panel the bands changed by the presence of potassium. 





for the optimized configuration using standard procedures in 
DFT. Fig. [3]gives some details of the Brillouin Zone form and 
orientation relative to the crystal lattice. Also the path used to 
calculate the band structure is shown. Fig. |4]gives results cor- 
responding to the LDA clean picene crystal (upper panels) and 
doped Kapicene in the LDA herringbone crystalline structure 
(lower panels). The most significant feature of picene elec- 
tronic structure is the opening of a wide gap (about 2.5 eV in 
our approximation) separating occupied (valence) states from 
empty (conductance) states. The upper right panel shows the 
bands that are close to the chemical potential. This result is 
in accordance with the insulating character of picene although 
the gap is underestimated in the local density approximation 
(LDA). After adding six potassium atoms per cell and assum- 
ing that the theoretically proposed crystalline structure is cor- 
rect, the material presents the electronic structure shown in the 
lower panels of Fig. [4] Both valence and conduction bands 
are globally similar to the previous ones except for two new 
peaks appearing in the valence band, the first just above -30 
eV and the second around -13 eV. They correspond to the fully 




FIG. 6: (Color online) Diffusion barrier (eV) for one K atom dif- 
fusing inside the frozen clean picene unit cell. Configurations for 
the initial and final states are shown in the upper and lower panels, 
respectively (hydrogens have been removed from the picture for the 
sake of clarity). 
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FIG. 7: (Color online) Diffusion barriers (eV) for two K atoms mov- 
ing in the same direction (triangles pointing down, continuos lines) 
and in the opposite direction (triangles pointing up, dashed line). 
Configurations for the initial and final states are shown in the upper 
and lower panels, respectively. 



occupied 3s and 3p states coming from potassium atoms. The 
valence electrons of K atoms have been transferred to the con- 
duction band as the Fermi level position indicates (we have 
checked that the integral below the occupied part of the con- 
duction band is precisely six). Consequently, the potassium 
4s states lie well above in the conduction band (but are not 
shown in the Figure) and are completely empty. A closer look 
to the states that are occupied by K electrons in doped picene 
is given by Fig. [5] It is clear that the dispersion of the bands is 
larger in the doped material and consequently, the separation 
to the valence bands is smaller. Notice also that dispersion 
along FQ and ZT directions (both parallel to k3 axis (c* axis) 
is more or less similar. Although the assignment of wave- 
functions weight to particular atoms or molecules is beyond 
the scope of this contribution, we believe that extended states 
of the doped material are similar to pristine picene wavefunc- 
tions but their energies are shifted by the new electrostatic po- 
tential induced from the transference of electrons from potas- 
sium to organic molecules. Therefore, assuming that 4s K or- 
bitals do not have a significant weight in the conduction band 
states, the existence of several weak CH/tt bonds 18 19 among 
perpendicular picene molecules can explain the existence of 
extended states. On the other hand, the smaller H-H overlaps 
should provide the smaller derealization along the c axis. 

Once changes in the electronic structure of picene due to 
the doping by K have been analysed, some conclusions can 
be inferred. First, the introduction of potassium gives rise 
to an ionic system in which long range Coulomb forces con- 
tribute an important percentage of the whole cohesive energy. 
Second, the broadening of the molecular levels of individual 
picene molecules likely comes from the original edge to face 
interactions of any herringbone arrangement plus the smaller 
interlayer coupling mediated by H-H overlap. Consequently, 
the doped system presents extended states covering the whole 
crystal, but showing a reduced dispersion along the ]<3 axis. 
This characteristic implies a smaller electric conductivity par- 
allel to the c axis in comparison with conductivities in the a 
and b directions. Measurements mentioned by Mitsuhashi et 
al. (Ref. f° r doped pentacene share this property although 
physical grounds are quite different. In the last case, conduc- 
tivity shows a two-dimensional character because doping oc- 
curs in a-b planes allowing extended two-dimensional wave- 
functions for the alkali metals. 



C. Intercalation of K 

Intercalation of K atoms inside the picene crystal is crucial 
to modify some properties of interest, e.g. superconductivity. 
In this subsection, we use our structural and electronic calcu- 
lations to compute diffusion barriers for K atoms moving in- 
side the unit cell. The starting point are local equilibrium ad- 
sorption sites computed for a common unit cell that we choose 
to be the experimental one for clean picene 2 ; this is a reason- 
able choice because of the low amount of potassium involved 
in these calculations. Barriers are estimated from a Linear 
Synchronous Transit (LST) transition state search, followed 
by a Quadratic Synchronous Transit (QST) refinemenPil 
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TABLE I: Barriers from the initial (Bi) and final (Bf) states, and 
energy of reaction (AE = Ef — Ei) for one and two K atoms dif- 
fusing in the unit cell (primarily along the c direction). For two K 
atoms two different cases have been considered: a symmetric (anti- 
symmetric) diffusion mode (S or A, respectively) where both atoms 
move along the same (opposite) directions. 

""System I B/ B F AE 
OneK 0.13 0.02 0.11 

TwoK(S) 0.61 0.33 0.28 

Two K (A) 0.47 0.44 0.04 

The global energy minimum for a single K in the picene 
unit cell is close to the central part of the channel defined 
by four picene molecules (Figs. [1J and [6]). We compute the 
diffusion barrier to move the K atom along the channel (the 
displacement introduces new forces making the atom to get 
closer to one of the picene molecules by 0.6 A). The energy 
of reaction when moving between these two local adsorption 
sites is 0.1 1 eV and the highest barrier (from the initial posi- 
tion to the transition state) is 0.13 eV (Table [I]). This barrier 
is small enough as to predict a large diffusivity of K inside 
the picene crystal: we estimate this value from We~ B / kBT , 
where W is a typical attempt frequency that we take as the 
frequency of a plausible phonon (e.g., 1 THz), B is the diffu- 
sion barrier for this particular hop from site I to site F, T is 
temperature and ks is the Boltzmann constant. This simple 
estimate shows that the potassium atom is making about 10 9 
hops per second at room temperature, making macroscopical 
diffusion in a reasonable amount of time possible. 

For two K simultaneously diffusing in the unit cell many 
more configurations need to be explored making an extremely 
complex energy landscape. We show in Fig. [7] two repre- 
sentative cases: starting out from the global minimum config- 
uration where both K atoms stay near the center of the unit 
cell, both atoms diffuse either along the same or opposite di- 
rections (mostly, both along the +c, or the +c and — c axis 
respectively). The lower barrier corresponds to the K atoms 
moving away from each other: 0.47 and 0.44 eV from the ini- 
tial and final states respectively. For the concerted diffusion, 
barriers are 0.61 and 0.33 eV respectively, and the initial and 
final states have a larger energy difference: 0.28 eV (Table|I|. 



IV. CONCLUDING REMARKS 

Systematic use of DFT allows the prediction of a plausi- 
ble crystalline structure for heavily doped picene (K^picene). 
Potassium is absorbed in the square prism space defined by 
four equivalent picene molecules. The original herringbone 
structure formed by the organic molecules is only weakly de- 
formed. This justifies the similarity of the X-ray powder spec- 
tra of the doped material and the pristine picene. Nevertheless, 
our predicted structure do not quantitatively reproduce the ex- 
perimental spectrum, mainly because the larger angle peaks 
originating in K are not present in experiment. The metal- 
lic character of the doped material comes from the occupa- 
tion of picene wavefunctions due to the transfer of electrons 
from potassium to the organic molecules. The presence of K 
cations induce a modified electrostatic potential that broad- 
ens the lower part of the conduction band without apprecia- 
bly changing the wavefunctions based on molecular orbitals 
of picene. Therefore the reason for the intralayer dispersion 
of the bands is the CH/ ir interaction responsible for the her- 
ringbone PAH arrangement. We should comment that this de- 
localization effect will be lost for separated parallel picene 
molecules, providing an important argument against alterna- 
tive laminar structures. In view of these results, we assume 
that superconductivity might be triggered by electron-phonon 
coupling within individual picene molecules. 
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